Monte Carlo simulation method for Laughlin-like states in a disk geometry 
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We discuss an alternative accurate Monte Carlo method to calculate the groundstate energy and 
related quantities for Laughlin states of the fractional quantum Hall effect in a disk geometry. This 
alternative approach allows us to obtain accurate bulk regime (thermodynamic limit) values for 
various quantities from Monte Carlo simulations with a small number of particles (much smaller 
than that needed with standard Monte Carlo approaches). 

PACS numbers; 73.43.-f, 05.30.Fk, 71.70.Di 



INTRODUCTION 



The discovery of the fractional quantum Hall effect (FQHE) has stimulated extensive studies on the properties 
of two-dimensional (2D) quantum many-electron systems in a strong magnetic field [ij. It is now understood that 
the FQHE represents the condensation of nearly 2D electrons subject to a strong perpendicular magnetic field (at 
low enough temperatures and low enough amount of disorder) into an incompressible quantum fluid formed at some 
specific filling factors. A neutralizing positive charge density is present to preserve overall charge neutrality and, 
to lowest approximation, can be thought of as a uniform positive density in the same plane as the 2D electrons (in 
reality the positive charges are the ionized donors that are roughly distributed randomly a spacer thickness away and 
produce a small amount of disorder and an overall constant shift in the energies). 

The strong magnetic field quantizes the electrons 's motion on the plane and quenches the kinetic energy of each 
electron to a discrete set of Landau levels (LL) separated by the relatively large cyclotron energy hujc — TieB/m*, 
where — e(e > 0) is the electron charge and m* is the effective mass of electrons in the semiconductor (to* ~ O.OTtoe in 
GaAs, where me is the bare electron's mass). In addition, the Zeeman splitting spin-polarizes the electrons rendering 
them effectively spinless. In each LL, p^L — l/27r^g (^o = ^/fi/eB is the magnetic length) electrons per unit area can 
be accomodated. It is evident that for large enough magnetic fields only the lowest LL (LLL) is occupied (and only 
the lowest spin sub-band), and if the occupation is not complete the system is highly degenerate. At particular filling 
factors V = p/ plLi electron interactions lead to highly correlated states which exhibit an excitation gap and result in 
the observed FQHE. These particular filling factors form a hierarchy, of which the simplest FQHE states have filling 
factors V — 1/to with to = 3 and 5. 

For filling factors of the form v = Ijin (n\ odd) the unnormalized Laughlin [5] trial wave function for N electrons 
can be written as: 

^^(Zi ...Zn)= X{{Zr ~ Z,T n ' (1) 

z<j i— 1 

where Zj = Xj + iyj is the position of the j-th electron in complex coordinates. This wave function gives an excellent 
description of the true ground state of the electrons for to = 3 and 5. For to > 7 the electrons tend to form a Wigner 
crystal |3] consistent with the experimental observation [4] that the FQHE does not occur for filling factors i' < 1/7. 
Since the Laughlin wave function lies entirely in the LLL, the expectation value of the kinetic energy per electron 

^ = h.^ , (2) 

is constant and becomes irrelevant, therefore the only important contribution of the quantum mechanical Hamiltonian 
H = K + V originates from the total potential energy operator: 



V^Vee + Veb + Hb , (3) 



where 

TV 9 N 



i<j 



2 



are the electron-electron, electron-background and background-background potential energy operators respectively. 
Here we have assumed a simple geometry approprite for the circular symmetry of the Laughlin wave function, where 
a positive background density po = ^/{'^t^^q) is spread over a disk fi^r of radius Rfq = Iq\J2N/v (i.e. it cancels the 
electronic density in the thermodynamic limit and makes the system neutral for all N). 

Numerous techniques have been employed to calculate the expectation value of the potential energy per particle 
{V)/N [see Eq. ^] in the Laughlin state [Eq. ((T|)]. For example, Laughhn 0] initially employed the hypernetted-chain 
method to estimate the value of this correlation energy with a ~ 1% accuracy); and various standard Monte Carlo 
(MC) schemes have been proposed, all of which are essentially exact in the thermodynamic limit. 

An excellent description of a standard MC computation of the potential energy and other relevant quantities in a 
disk geometry is given by Morf and Halperin [H| . Spherical geometries are also used quite often since the convergence 
to the thermodynamic limit is quicker because boundary effects are eliminated 

Although considerable more care is needed in the disk geometry to eliminate boundary effects when extrapolating 
(necessarily) finite- iV results to the thermodynamic limit (in particular due to the long-range nature of the Coulomb 
potential), there are cases in which the spherical geometry is either inconvenient, or plainly incompatible with the state 
under consideration (for example for the study of possible quantum Hall nematic phases 0] for which considerable 
topological defects would be generated at the poles of the sphere). 

Furthermore, the value of the correlation energy in the thermodynamic limit is not easily extracted from standard 
MC simulation (see Sec. ) data, since the limit is approached very slowly, with corrections of the order of 1/Vn, 
requiring repeated calculations for various N and a careful extrapolation of the results to the N —^ oo limit. It is, 
therefore, highly desirable to explore methods that would expedite the extrapolation to the thermodynamic limit. 

In Sec. we describe, for the sake of completeness, the procedure used in the standard MC approach. Section 
describes an alternative method that converges to the thermodynamic limit considerably faster. We discuss our 
results in Sec. . 



THE STANDARD MONTE CARLO APPROACH 



In the standard MC approach one considers the calculation of the expectation value of the potential energy oper- 
ators as given in Eq. The background-background interaction potential poses no problem, it can be calculated 
analytically and is found to be 



In order to compute the expectation value of the electron-background interaction potential one conveniently writes 
it as 

Veb = y]vebiri) ; VebiYi) = -Po / d^r- -, (6) 

where Veti'^^i) is the interaction potential of a single electron at position fi with the uniform positive background in the 
finite disk. Such electron-background interaction potential depends on the ratio ri/R^, where — \ri\ is electron's 
distance from the center of the disk and Rn is the radius of the finite disk and can be expressed as 



, r , y)Ji{y) 



— , a; < 1 



vM = -V2^ Fin/RM) - ; Fix)^ dy ^^^^ ^ <^ , (7) 

'o Jo y \ ^p,(i 1-2- 1 



2x 



x> I 



where Jnix) are Bessel functions of order n, E{x) is the complete elliptic integral, and 2Fi{a,b; c; z) is the hyper- 
geometric function. Figure [T] plots the function F{x). It is interesting to note that F{0) = 1, F{1) = 2/tt and 
F{x) ^ l/{2x) for X ^ 1, as expected for the Coulomb potential far from a charged disk. Although F{x) can be 
expressed analytically, it generally preferable to store it in a table, and interpolate it for all x for all the calculations 
that follow. 

While most electrons stay within the confines of the neutralizing background (i.e. x < 1), electrons near the edge 
may spread outside the disk to some extent (although it is extremely unlikely that they will spread to more than a 
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FIG. 1: The electron-background function F[x). Here x — t/Rn, where r is the distance of an electron from the center of 
the disk and _Rjv is the radius of the finite disk filled with neutralizing positive background. The dotted line shows the l/(2x) 
asymptotic dependence. Most electrons sit in the x <1 region. 



few magnetic lengths from the edge). The expectation value of the electron-background interaction potential during 
the MC simulation can then be calculated using 



{Vel 



N 




(8) 



Finally, the expectation value of Vee is accordingly given by 



N 




(9) 



In the usual Metropolis MC method [s^ , the expectation value of an operator can be computed by averaging the 
value of the operator over numerous configurations {ri, . . . , r/v} of the many-body system that obey detailed balance, 
that is, the probability ratios between pairs of discrete configurations are related by the ratios of the probability 
distribution for the system [in this case |4',„(zi, . . . , zjv)^, see Eq. ([T])]. Usually several million configurations are 
used for each N and the results are extrapolated to the thermodynamic limit by considering a sequence of various 
increasing iV-s. 

A MC step (MCS) consists of attempts to move one by one all the electrons of the system by a small distance of 
order A in a random direction. After each attempt (to move the i-th electron from 7^°'^ to r'j"'^'", the probability 
ratio between the "new" state and the "old" state is then computed: 



„old 



■riv)| = 



cxp 



N 



new 1 2 



In |r, 



„oId|2\ 



exp 
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2P 



(10) 



In the usual Metropolis scheme if this ratio is bigger than a uniformly distributed number in the [0,1] range the 
attempt is accepted, otherwise it is rejected. The parameter A is adjusted so that the acceptance ratio is close to 
50%. After attempting to move all electrons (one MCS), the electron configurations are then used to calculate the 
operator under consideration. Averaging over numerous MCS-s converges gradually (as l/\/number of MCS) to the 
desired expectation value. Normally it is convenient to disregard numerous (several thousand) initial configurations 
to reach a good "thermalization" before the averaging begins, which significantly reduces the expurious effects of the 
somewhat arbitrary initial configurations. All the results that we report here were obtained after discarding 100,000 
"thermalization" MCS-s and using 2 x 10^ MCS-s for averaging purposes. 
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TABLE I; Correlation energy per particle in the Laughlin state for filling factors u — 1/3 and 1/5. These results were obtained 
after a standard Monte Carlo simulation in a disk geometry. Energies are in units of e^/lo- 



N 


m=3 


m=5 


4 


-0.38884 


-0.32159 


16 


-0.39766 


-0.32328 


36 


-0.40129 


-0.32446 


64 


-0.40323 


-0.32510 


100 


-0.40445 


-0.32550 


144 


-0.40521 


-0.32577 


196 


-0.40579 


-0.32594 


400 


-0.40675 


-0.32624 



In Table U we show the correlation energy per particle for finite systems of N electrons in the Laughlin states m — 3 
and m = 5 obtained using the standard MC method described above. The results are rounded in the last digit. 

To get the the bulk (thermodynamic estimate) of the correlation energy per particle one needs to perform a careful 
extrapolation of the results as illustrated in Fig. [2] where we show the correlation energy per particle for states v = 1/3 
and 1/5 plotted as a function of 1/^/N for systems with N = 36, 64, 100, 144, 196 and 400 electrons. 

We fitted the energies of Table H] for iV = 4, 16, 36, 64, 100, 144, 196 and 400 electrons to a polynomial function as 
reported in Ref. \^ and obtained: 



These interpolation lines are used to estimate the correlation energy per particle in the thermodynamic limit (the 
first term in each of the parentheses). Our results for thermodynamic limit are similar to those found in Ref.jol, 
-0.4100 ± 0.0001 and -0.3277 ± 0.0002 (in units of e'^/lo) derived with the use of the pair correlation function 
evaluated from MC simulations with up to = 256 electrons and generating as many as 5 million MC configuratons. 
Note how slow the convergence is, although the extrapolation to ^ oo seems unambiguous it is still time consuming 
and cumbersome ( even for N — 400 the error is still circa ~ 1 % ). 



The one-particle density 

Other physical quantities of interest that may be readily computed are the single-particle density function and the 
pair distribution function. Given that the Laughlin wave function describes an isotropic liquid and is rotationally 
invariant, the single-particle density depends only on the radial distance from the center of the disk. We may compute 
the single-particle density by counting the number of electrons iV;(A7') found in several 2D shells of width Ar centered 
around a discrete set of distances to the center ri = (I + ^)Ar (^ = 0, 1, . . .): 

where r2/(Ar) = 7r(Ar)^[(Z + 1)^ — P] is the area of each 2D shell. In the Ar the quantity computed corresponds 
unequivocally to the electron density 

p{r) ^ (j2S{r - (14) 

The computation of the single-particle density in the Laughlin state, indicates a significant nonuniformity near the 
boundary (see Fig. [3]). As the number of electrons increase, a significant portion of the system becomes uniform 



5 




-0.322 - 



A 
> 

V 



-0.324 



-0.326 



-0.328 - 



-0.33 I ' ' ' 1 

0.05 0.1 0.15 0.2 

1/sqrt(N) 

FIG. 2: Monte Carlo results in disk geometry for the Laughlin state at i/ — 1/3 (top panel) and u = 1/5 (bottom panel). 
The potential energy per particle, {V)/N is plotted as a function of for systems with A'^ = 36, 64, 100, 144, 196 and 400 

electrons. Full circles: correlation energies calculated by the standard method described above, the full line is a least-square fit 
[Eq. (|lip and (|12p] used to extrapolate to the thermodynamic limit. Energies are in units of e^/lo- 



as expected. Note, however, that the non-uniformity near the edge always persists. This behavior can be used to 
characterize which electrons are "in the bulk." 




5 10 15 20 25 30 35 40 

FIG. 3: One-body density function, p{r)/po, in the Laughlin state v = 1/3 as a function of the distance r/lo from the center 
of the disk for systems with A'^ = 64, 100, 144 and 196 electrons. Note the persistence of an "edge region" of finite width and 
the development of a "bulk region" for large A''. A discretization interval Ar — 0.05 lo was used. 



The pair distribution function 

Another important quantity related with the trial wave function is the pair distribution function, which corresponds 
to the conditional probability density to find an electron at a distance r from another electron. For any homogeneous 
and isotropic liquid with uniform density po it is defined as 



/ AT JV \ 
\i=l I 

Following the same procedure as above, we discretize in concentric shells around the i-th electron and count the 
number of electrons Ni{l^r) in each shell, which should give g{r) as Ar — + according to the following equation: 

1 1 / ^ ^ \ 
f^^^^^^) = NnK^){T.Y.^^i^r)) . (16) 



= 1 j^i 



It is evident that electrons near the edges of the system could contribute expuriously to these sums as their "surround- 
ings" are considerably different than those at the bulk. To eliminate as much as possible any boundary effects, it is 
convenient Q to consider only (for the "z-electrons" above) the electrons that are within a small circle of radius i?i 
around the origin. If iVi is the average number of electrons that are within this small circle, then the approximation 
g{ri) for the pair distribution g{ri) is: 

I N^ N \ 

\4 = 1 j^t I 



where in this expression one is considering the pairs between any electron i (i = 1, . . . , Ni) lying inside the circle of 
radius with all other electrons j {j — 1, . . . ,i — l,i + 1, . . . , N) that may lie either inside, or outside that circle. 
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This guarantees that the evaluation of g{ri) involves only pairs, where at least one member lies inside a circle of radius 
Ri around the origin, where correlations are believed to be close to those in the bulk of an infinite system. 

Fig. m show plots of the pair distribution function for the states v — 1/3 and 1/5 for systems with N = 
4,16,36,64,100,144 and 196 electrons. For our MC simulations we chose Ri = 0.25 R^ and a discretization in- 
terval Ar = 0.05 If). Note the gradual decay at large r which reflects the finite size of the system. 




FIG. 4: Pair distribution function for f = 1/3 and 1/5 obtained by a standard Monte Carlo simulation in disk geometry for 
systems of TV = 4, 16, 36, 64, 100, 144 and 196 electrons. 



The determination of the pair correlation function for a given finite TV tends to be quite time-consuming but provides 
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for an alternative way to compute the correlation energy per particle in the thermodynamic limit Q by using the 
formula: 

which is valid in the limit of an infinite system. Although the pair distribution function is obtained from a system 
with a finite number of particles, one can calculate the thermodynamic value of the correlation energy per particle to 
a very good accuracy by using the slightly modified formula 



N 2 

in connjuction with the normalization condition 



d'r - [g{r) - 1] , (19) 



po / d^r [g{r) - 1] = -1 , (20) 

which defines an upper cuttoff Rcut ■ This approach produces good estimates for the thermodynamic correlation energy 
per particle as long as g{r) is able to reach its asymptotic value [g{r) ~ 1]. Reasonable results can be achieved even 
for systems of > 36 electrons. 



AN ALTERNATIVE METHOD 



In the standard MC approach one needs to calculate the expectation values of various quantities several times for 
different N in order to extract the thermodynamic estimate from the data by performing a fit (and taking 

the limit 1/VN ^ , as shown, e.g. in Fig. [51 It is highly desirable to obtain estimates for the bulk regime (in 
thermodynamic limit) without needing to perform the above analysis. In the following we describe a method that 
allows us to obtain results consistent with the bulk regime, by doing simulations with only a finite (relatively small) 
number of particles. The method although approximate, yields very accurate results corresponding to the bulk regime 
even when simulations are performed for a small number of electrons. The estimates are very stable over a wide range 
of N and the technical application of the simulation is less involved. 

In order to obtain reliable estimates for the bulk regime, we need to exclude from consideration the boundary- 
affected outer region of the finite disk. In a standard MC simulation the electrons are distributed all over the 2D 
space, at a given instant it is obvious that the electrons close to the central region of the disk resemble to the bulk 
regime better than those close to the boundary. However during the simulation each of the previously "bulk" electrons 
moves around the whole disk therefore the correlation energy of such electron with the other electrons is not a good 
estimate of the correlation energy in the bulk regime. 

The core of the method proposed here is to consider an electron pinned to the center of the disk which, by 
construction is the point which most closely resembles the bulk of the system. Therefore, if we are able to derive 
results where only the correlation energy of that particular electron with the rest (away from edge) is involved, we 
anticipate that such estimates should approximate the bulk regime much more accurately than other methods, and 
as a result thermodynamic limit values can be achieved even in a system with a relatively small number of electrons. 
One has to recall that for a finite system, the Laughlin wave function describes an incompressible system of strongly 
correlated ellectrons with uniform density po only in the bulk region (central part of the disk not very close to the 
boundary), while close to the boundary (where the density of electrons falls to zero) the fluid becomes compressible 
and there is a deviation of the electron density from its constant value in the bulk. 

In our MC simulations, we consider a Laughlin-like state [Eq. ([T])] in which one electron is pinned at the position 
zo (we consider zq = 0), and N' = N — 1 electrons are free to move (i.e. in a typical MC step): 

^'^{zo,zi,...ZN') ^'i'„-,izi,...ZN') Y[izj-zo)"'e . (21) 

j=i 

When an attempt is made to move electron i (i = 1, . . . , N') from r°''^ to rf™, the MC probability ratio is given by 

|^'(ro,...r— ...r^OP _ |vl/(ri, . . . rf- ■ ■ ■ r^QP ™(i„ |r --ro^-m ir-^-roP) (^o) 
|vI/'(ro,...r°id...r^,)|2 |vl/(ri, . . . rf ^ . . . r^,)|2 ■ ^ ^ 
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TABLE II: Correlation energy per particle in the Laughlin state for filling factors v = and 1/5, obtained via a Monte Carlo 
simulation in a disk geometry using the method of pinning one electron at the center of the disk. Energies are in units of /Iq. 



N 


m=3 


m=5 


4 


-0.38187 


-0.30157 


16 


-0.40898 


-0.32722 


36 


-0.40895 


-0.32637 


64 


-0.40909 


-0.32665 


100 


-0.40955 


-0.32738 


144 


-0.40936 


-0.32732 


196 


-0.40953 


-0.32734 


400 


-0.40954 


-0.32735 



Since the electrons are identical we need only consider the average correlation energy of one specified electron with 
the rest of the electrons to compute the electron-electron interaction energy, therefore to this level of simplification 
the electron-electron energy per particle is obtained by considering only the interaction between the pinned electron 
and the rest of Ni other electrons contained within a inner disk of radius Ri < Rn where the electron density is 
approximately equal to the bulk value pQ. We found that a reasonable choice for the radius of inner disk that excludes 
the edge electrons is Ri = 0.75 Rn, therefore this value was adopted in all the following simulations. 

To be consistent with the above procedure also the disk's electron-background and background-background energy 
should be calculated within the same degree of simplification (see Appendix ). It is, therefore, useful to first calculate 
the one-body density function (see e.g. Fig. [3]) in order to determine an optimal Ri for future use. If Ri is reasonably 
large, we expect the correlation energy calculated in this way to closely correspond to the desired correlation energy 
per particle in the thermodynamic limit (see Appendix): 

iV - 2 \^ \n-ro\/ \ 2m lo ' ^ ' 

As in the standard MC method case, our MC runs consist of 100,000 discarded equalibration MCS-s followed by 
2 X 10^ MCS-s for averaging purposes. 

In Table |TT] we show the correlation energy per particle for finite systems of N electrons and Laughlin states m = 3 
and m — 5 calculated by pinning an electron at zq = as described above [see Eq. ([22]) ]. and using Eq. ((23)) . The 
results are rounded in the last digit. 

In Fig. [5] we show the potential energy per particle for v = 1/3 and 1/5 computed from the alternative method and 
plotted as a function of 1/^/N. For the sake of comparison we also plot the results from the standard method (Sec. ). 

It is striking to note how much faster the alternative method converges to the thermodynamic limit. Differently 
from the standard MC approach, the alternative method that we introduced does not need to have the data points 
least-square fitted to get the thermodynamic limit value. One merely needs to choose a big enough (for instance 
N — 196) and do a MC run that typically will generate bulk (thermodynamic limit) results to an excellent degree of 
accuracy. 



The pair distribution function 

Using the same ideas presented above (keep one electron pinned at zq — 0) the pair distributiuon function is very 
easily calculated. The essence of the method consists in measuring the one-particle density excluding the pinned 
electron which is obviously the pair distribution function instead of considering all the possible pairings between the 
electrons. As before, the advantage of this method, besides its simplicity, resides in the fact that this electron, being 
the farthest from the edge, is in an environment closest to that in a bulk system. Therefore, one has to calculate: 
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FIG. 5: Monte Carlo results in disk geometry for the Laughlin state at = 1/3 (top panel)and 1/5 (bottom panel). The 
potential energy per particle, {V)/N is plotted as a function of for systems with A'^ — 36,64,100,144,196 and 400 

electrons. Full circles: correlation energies calculated by the standard method (Sec. ) the full line is a least-square fit [Eq. (|lip 
and (|12[) ] used to extrapolate to the thermodynamic limit. Empty circles: correlation energies calculated by the alternative 
method described in Sec. , the dashed line is a visual aid indicating the thermodynamic limit. Note how the thermodynamic 
limit is approached faster in the second method. Energies are in units of /Iq. 



where Ni{Ar) is the number of electrons found in the 2D shell f2i(Ar) with distance range (r;, r; + Ar) from the pinned 
electron at the center of the disk. In Fig. [6] we show plots of the pair distribution function for the states 1^ — 1/3 and 
= 1/5 for systems with N = 4,16,36,64,100,144 and 196 electrons obtained with a choice of Ri — 0.75 Rn and 
using a discretization interval Ar = 0.05 Zq- 
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FIG. 6: Pair distribution function for the state v = 1/3 and 1/5 obtained after a Monte Carlo simulation in the disk geometry 
for systems of iV = 4, 16, 36, 64, 100, 144 and 196 electrons with one electron held fixed at the center of the disk . 



SUMMARY AND DISCUSSIONS 



We implemented an alternative Monte Carlo method to calculate the properties of Laughlin states of the fractional 
quantum Hall effect in the thermodynamic limit while using a very small number of electrons. The key point of this 
method is the pinning of an electron in the center of the disk, so that the potential energy and correlation functions 
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calculated through the pinned electron accurately represent the bulk (thermodynamic limit) even in a relatively small 
system. The idea is quite general and, in principle, can be applied to any system as far as the main concern is the 
calculation of correlation effects such as t he p otential energy, the pair distribution function, etc. For systems such 
as 2D electronic one-component-plasmas [10], composite fermion states described by the Jain's unprojected wave 



function [ll[ etc, that is all that matters. Obviously such alternative method can always be used to calculate the 
potential energy and related quantities of other more diverse systems, with the the kinetic energy calculated in the 
standard way whenever applicable. By using this alternative method we analyzed the properties of the Laughlin states 
corresponding to filling factor v — l/i and 1/5 by performing Monte Carlo simulations in disk geometry for systems 
with up to = 400 electrons. The correlation energy per particle and the pair distribution function computed in this 
approach are compared to corresponding bulk-regime values obtained via a standard Monte Carlo simulation in disk 
geometry, where a careful extrapolation in thermodynamic limit has been performed. We find that such approach 
allows us to obtain accurate bulk regime (thermodynamic limit) values for various quantities using a modest number 
of electrons ( even for iV = 16 the error is less than 0.1 % ). 
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APPENDIX 

Let us consider a system of N interacting 2D electrons coupled to the positive neutralizing background that fills 
a finite disk and guarantees overall charge neutrality. The expectation value of the electron-background potential 
energy per particle can be written as: 

(Kb) Po f ,2„ N /■ j2„ e 



/ <fr,p{v,) f d\ , (25) 



N N 

where p{vi) is the one-body (single) electron density function given by 

P(ri) - LJ ^ . / pin) = N . 26 

J d^ri . . .d^rAT |^'(ri . . .rAr)r J 

One notes that when the one-body electron density becomes uniform, p(ri) w po and the system is sufficiently large 
so that most of the electrons are to be found inside the finite disk then —{Veb)/N « 2{Vi,b)/N, therefore we would 
have: 

jVeb) + jVbb) ^ 1 (Veb) 

N - 2 N ■ ^ '> 

By making a preliminary calculation of the one-body density function we could estimate the radius, Ri < Rn of an 
inner disk where the electrons have supposedly uniform density therefore we could argue that all these electrons inside 
this inner disk are in the bulk regime. If Ni is the number of electrons within this reference circle of radius Ri in 
addition to the electron pinned at the center of the disk then a good estimate for the potential energy per particle in 
the thermodynamic limit is obtained from the quantity: 

where {Vhb)/{Ni -f 1) is the background-background energy per particle of a positive charge that exactly neutralizes 
the charge of Ni + 1 electrons. Since the electron-electron potential energy per particle was obtained to a level of 
simplification where only the interaction between the pinned electron and other iV^ electrons was considered, then the 
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whole potential energy per particle should be calculated in the same grounds too, namely, using as reference only the 
interaction energy of the pinned electron with the positive background: 

As a result the correlation energy per particle in the thermodynamic limit can be written as in Eq. (|23p. Note that 
use of relation ~(Veb)/N « 2{Vbb)/N to express + {Vbb))/N in Eq. ip?)) in terms of {Vbb)/N is inaccurate in view 

of the approach adopted in our method therefore should be avoided. 
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